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A localized wave packet exhibits under certain conditions zitterbewegung. A similar phenomenon 
occurs for the chirality of a massive particle. In the case of massive neutrinos, since they are detected 
via the V-A weak charged currents, this oscillation may even explain the "missing" solar neutrino 
experiments. The neutrino remains a mass eigenstate but contains an almost sterile right-handed 
component. This qualitative discussion opens up a number of interesting physical considerations. 
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We are interested in chiral oscillations as a possible explanation of the missing solar neutrino problem because, 
neutrinos with positive chirality are decoupled from the neutrino absorbing charged weak currents. In recent years 
there has been enormous interest in the problem of the missing solar neutrinos . With the inclusion of the latest 
experimental results, approximately a third of the expected solar electron neutrino flux, within the energy ranges 
' measured, is missing in earth bound experiments The related search for neutrino masses was stimulated by the 

lO . mass oscillations invoked as the most promising hypothesis to explain the missing neutrino data. If the neutrinos have 
mass then it is possible to postulate mixing between known neutrino species (electron, muon and tau) Q or between 
the known neutrinos and higher mass neutrinos including sterile neutrinos ||. Such assumptions together with the 
violation of individual (e, fi, r) leptonic number permit oscillations between neutrino species. This hypothesis is 
extensively studied in the literature , and is characterised by the parameterization with the differences between the 
squared values of neutrino masses and, of course, by the mixing angles. As we shall see below, our chiral oscillations 
""^ ' are quite different, and do not involve any mixing angles. 

Let us forget the neutrino for a moment and recall some basic facts about zitterbewegung . In first quantization, 
the relativistic Dirac equation predicts the existence for the spin i particle, under certain conditions, of a time 
Qh! dependence (oscillation) in the average of the position variable (the velocity eigenvalues associated with this oscillation 
^ are ±c even for a massive particle). This phenomenon, known as zitterbewegung, is a consequence of the non 
commutation of x with the Hamiltonian H . Since our calculation for chirality follows an almost identical line of 
^ i reasoning, we shall not give here any details of the demonstration of zitterbewegung. We recall the possibility to 
'• av °id zitterbewegung by redefining the space variable via the Foldy-Wouthuysen transformation However, this 
introduces a very complicated position variable and is in contrast with the justification of the Darwin term in atomic 
physics [ p7o[ | . In fact, the existence of zitterbewegung is no stranger than the existence of negative energy solutions, 
on which it depends. We wish to recall that zitterbewegung is only manifest for wave functions with significant 
interference between positive and negative frequencies. In the following we shall talk of frequencies in plane waves 
and not of energies. Frequencies and energies coincide only for plane waves . For bound states and localized production 
the Hamiltonian is never the free particle one and frequency and energy no longer coincide. For example, the energy 
eigenstates of the hydrogen atom are not plane wave solutions and indeed contain both types of frequency, but, non 
the less, they represent definite energies. 

For an electron, model calculations pd| ] predict an oscillating zitterbewegung frequency of around 10 21 cycles per 
second, far to rapid at present for any hope of direct measurement. However, we mention as an aside that if the 
neutrinos are massive, one of them, probably the electron neutrino, is the least massive of all known fermions and as 
a consequence its zitterbewegung frequency may even be within the capacity of direct experimental time resolution. 
Even if we cannot envisage any practical experimental test, given the neutrino's ephemeral properties. 
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The starting point for what follows is the observation that the creation of a particle is always a localized operation, 
notwithstanding the use in perturbation theory of plane waves (the asymptotic as well as unperturbed wave functions). 
For example a particle produced in the laboratory is obviously localized, sometimes its origin is measured to within 
a few microns in vertex detectors. Normally this fact is ignored or considered irrelevant. Nevertheless, the creation 
of a particle involves, even if for only short times, a non free Hamiltonian which fortunately we do not need to know 
for our present calculations. We shall simply assume the form of the wave function at creation as a gaussian, to ease 
calculations, hoping that this assumption at least approximates the true wave function. The average chirality of the 
particle then exhibits a time dependence similar to zitterbewegung, and it is this type of oscillation which we consider 
in this paper. 



II. CHIRALITY TIME DEPENDENCE 



If a spin ~ particle is massless and is produced in a chiral (7 s ) eigenstate, then its chirality is a constant of motion in 
addition to being Lorentz invariant. These properties no longer hold for Dirac particles with mass. The free particle 
Hamiltonian, 

H = -ia. ■ d + mj3 , (2.1) 

which we shall assume to represent the time- evolution operator for times subsequent to the creation of the spinor 
(t = 0), does not commute with the chiral operator 7 s because of the /3(= 7 ) term in H. Specifically, 

-zd t (7 5 >H[ff,7 5 ]}-2m<7° 7 5 > (2.2) 

where m is the mass of the spinor, assumed a mass eigenstate. The mean value on the right hand side of the above 
equation is given explicitly by 

(7°7 5 ) = J d 3 x ¥(a;)7 5 *(x) , (2.3) 
where $ may be expanded in terms of plane wave solutions in the standard way: 

a a (k)u a (k)e- ik * + b* a (k)v a (k)e ik * } , (2.4) 

J ^ W ' a=l 

with E which denotes the positive quantity 



E = k° = Vfe 2 + m 2 . 

The normalization of the spinors, 

v i( k ) v fj(k) = u+(k)up(k) = — S af3 

-V a (k)Vp(k) = U a {k)up{k) = S a /3 , 

and the normalization condition for ^f(x), 



impose the following constraint on the coefficients a a (k) and b a (k) 

d 3 k 



(2tt) 3 E 



E [M*)l a + IM*)l a ] =i- (2-5) 



Now, as for zitterbewegung, the formula ( |2.2| ) does not necessarily imply a physical effect. In it is easy to see that 
for wave packets made up of frequencies of only one given sign, the Dirac equation for the spinors u(k) and v(k) 

(ft - to) u(k) = u{k) (ft - m) = , 
()t + m)v(k) =v(k) (ft + m) = , 
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leads to a zero a vera ge in Eq. (|2.2|), To demonstrate this point, let us superpose only positive frequency plane waves. 



Computing Eq. (2.3), we find 

(7V) = (2^ ^ lJ ^ al{ q )a p {k)u a { q )^u p {k)e^- k > 

(2^)3 & <(fc)a/3(fcR(^)7 5 w/3(^) ■ (2.6) 



where E' = q° — \J q 2 + m 2 . Now as a consequence of the Gordon-like chiral identity, 

u a (q)"f 5 up(k) = — u a {q) (#7 5 + 7 5 #) u p {k) , (2.7) 

we have 

u a (k)j 5 ui3(k) = , 

and so we obtain a zero average for 7°7 5 . A similar proof works out for a superposition of only negative frequency 
plane waves. 

On the other hand if both frequencies are present the interference terms yield an oscillation in chirality, 

(7V) = j d " X (0 Wf ^E' 5 [ 6 «^) a ^ fc )^^)^^( fc ) e " %+ " )X ~ h ' C - 1 

]T [ b a (k)a p (k)U a (kh 5 u (k)e- 2iEt - h.c. ] , (2.8) 



(2tt) 3 E 2 

where k = (E, —ft). The Gordon-like chiral identity for cross terms, 

v a {q)^up{k) = v a (q) (^ 7 5 - 7 5 #) tt/J (fc) , (2.9) 

gives us, the non null equivalence 

v a (kh 5 u p (k) = -- v a Ckh°l 5 up(k) . (2.10) 
m 

From now on we will work with the following (Pauli-Dirac) representation of gamma matrices 
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In such a representation the spinor solutions are given by 



lie, (ft) 
Wo (ft) 

where conventionally 



77 i \ V 2 / 



2m / 1 CT 'fc' 



i /a / 2—^- -v/ 



2m 



(2.11) 



(2.12) 



1 \ / 



Xi = I , X2 - . 1 



Using Eqs. (2.11) and (2.12), we easily check that 



and so, combining this last result with Eq. ( 2.10 ), we obtain 



E E ~ ~ E 

v a (k)j 5 ut3(k) = v a (k)-/ j 5 Uf3(k) = — v a {k)vp{k) = 5 af} . (2.13) 

m mm 



Then Eq. (gj) becomes 

dtil 5 ) =J-0f E [ -^Eb a {k)a a {k)e^ Et + h.c. ] . (2.14) 



After performing the time integration, we find 

<7 5 >( t ) = (7 5 )(o) + / |J E t Uk)a a (k) (e^ Et l) + h.c. ] . (2.15) 

Even this formula does not necessarily imply a significant physical effect. For example, we shall demonstrate below 
that for a gaussian initial wave function and if the initial state has average chirality zero the oscillations of the ± 
chiralities cancel and there is again no overall oscillation. 

At this point we wish to discuss an apparent paradox. Now, for any plane wave solution of the Dirac equation with 
mass, the rest-frame wave function is always an equal mixture of both chiralities. This is easily seen by, 

i> = —y~ ^ + 2 i 3 

= ^++^-, (2-16) 
where "0+.- correspond to chirality ±1 respectively. Furthermore, we have that 

$i>± = +mip T , 
]iip± = -mip T , 

where +m appears for positive energy plane waves and — m for negative ones, and where k is the 4-momentum of the 
particle. Thus, in the rest-frame (fe = 0) we find, since (7 ) 2 = 1, that 

hM a = hM a • (2.17) 

Note that this result is not Lorentz invariant since a Lorentz boost is not a unitary transformation. Thus, while the 
cross section is Lorentz invariant the chiral probabilities are not. This seems to suggest that cross-section measurements 
are chiral independent. We seem to have an argument against the physical significance of chiral oscillations. 

The reply to this objection based upon the Lorentz invariance of the cross-section is simply that in any given 
Lorentz frame chiral oscillations are manifestly important because of the chiral projection form (V-A) of the charged 
weak currents. The chiral probability variations produced by Lorentz transformations (even if 7 s commutes with the 
Lorentz generators) are automatically compensated by the wave function normalization conditions and the Lorentz 
transformations of the intermediate vector bosons and other participating particles. 

This situation is similar to the apparent contradiction between the following facts. In all known processes the 
neutrino is produced as an almost left-handed helicity eigenstate because created dominantly in an ultra-relativistic 
state (where chirality ~ helicity) by the V-A interaction and hence with average chirality close to —1. But, in its 
rest frame the (plane wave) neutrino has average chirality zero. Viewed from the neutrino's rest frame there is no 
preference between + and — chiralities. Is this not in contradiction with the dominance of left-handed neutrinos in 
nature? No, because while there exists frames common to almost all neutrinos in which the chirality is almost pure 
— 1, there exist many, but non coincident frames for which +1 chirality dominates for each given neutrino. 

We conclude this Section with a discussion of another important question: in which frame shall we assume a given 
form of localization? A spherical gaussian assumption, suggested primarily by its simple integration properties, is 
certainly not frame independent. Thus if the sun is the relevant frame then we can assume, as we shall do below, that 
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the neutrinos are created as almost pure negative chiral eigenstates. Chiral oscillations then occur and this in turn 
implies that a measurement of the neutrino solar flux on earth will detect an "apparent" loss of neutrinos. 

Obviously, any neutrino created in the sun can reasonably be considered localized within the solar core at time zero. 
It is not even excluded that the relevant distance scale for localization may be much smaller such as the dimensions 
of the nucleons within which the decaying quarks current gives rise to the neutrinos. In this case however we would 
have to take into account the wave function of the parent hadrons. In the extreme hypothesis, that the distance of 
localization (at creation) of the neutrinos is the Compton wavelength of the intermediate vector boson W, we will see 
that the chiral oscillation is very small. 



III. A GAUSSIAN MODEL 



Let us assume that at creation the neutrino wave function is given (or approximated) by a gaussian function in x 
with range d for the corresponding probability distribution, 

V(0,x) = (ird 2 )~ 3/4 exp(-a; 2 /2d 2 ) w , (3.1) 

with the spinor w normalized by w + w = 1. Since we have already noted that both theoretically and experimentally 
the ultra- relativistic neutrinos are created in almost pure —1 chiral eigenstates (massive neutrinos cannot be rigorously 
in a chiral eigenstate), we assume for w the approximate form 

»-(£). ("> 

with ip a two component spinor, such that ip + ip = 1/2. 

Eq. (3.1) may be used to obtain the value of the coefficients a a (k) and b a (k), 

a a (k) = (47rd 2 ) 3 / 4 exp(-fc 2 d 2 /2) u+(k)w , 
b* a (k) = (4^ 2 ) 3 / 4 exp(-fc 2 d 2 /2) v+(k)w . 

Note that starting with a negative chiral eigenstate we find b* a {k) = —a a (k). Eq. ( 2.15| ) then yields the time dependent 
value of the average of 7 5 , 

(7%) = -1 + / |J E t <( fc K« (! " e " 2 ^) + h - c - 1 > ( 3 - 3 ) 
which after algebraic manipulations becomes 

<7 5 > (f ) = -l + NCO , (3.4) 
where NCO (Neutrino Chiral Oscillation) indicates the following integral 

NCO = / ^ (4W 2 ) 3 / 2 exp(-fc 2 d 2 ) [1 - cos(2£7i)] . (3.5) 

J (2w) i E z 

Alternatively we may quote the probability of finding the (1 — 7 5 )/2 component relevant for weak charged interactions, 
which coincides with its mean value: 

(-) = J d 3 x 9+(x) «f(x) . 

Explicitly, we have 



By changing the initial chirality to +1, 



H = 1 - \ NCO . (3.6) 
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we obtain b* a (k) — a a (k), and 



(+) = 1+ 2 NC0 • 

Consequently, we find the anticipated result that for average intial chirality zero, there is no oscillation. 
In spherical coordinates, Eq. (|3.5|) becomes, after angular integration, 



(3.8) 



NCO = 



A/2 



(mk) 2 

dk ^— f- exp(-kV) [1 - cos(2£i)l , 



(3.9) 



where k = |fc|. For future considerations, we prefer to rewrite Eq. (p.9[) after the following change of variable, p = kd, 



NCO 



rV2 



dp p 2 



2 d 2 



1 



eM-p 2 ) { 1 



2mt 



2 d 2 



1/2 



This equation is particularly simple when md may be considered large, 

NCO - 1 - cos(2mt) . 

Then the expectation value for the negative chirality average becomes 

. . 1 + cos(2mi) . , . 
(-) 7T " [ md » 1 ] • 



(3.10) 



The neutrino in this limit oscillates between negative and positive chiral eigenvalues. When in the chiral +1 cigenstate 
the neutrino only participates in elastic neutral current interactions and is thus in practice "missing" . To obtain a lower 
value for the neutrino mass we must determine the value of d. If our gaussian wave function and correspondent gaussian 
probability for neutrino production is identified with the luminosity [ ^2|JT3] | of the solar core, we find consistently a 
value of d of ~ 0.1 Rq (where Rq is the solar radius). This yields a neutrino mass > 3 x 10~ 15 eV. In any model 
of oscillations we must allow for the fact that the arrival times of the solar neutrinos varies. In a purely particle 
language (very short range localization) we would say that this depended upon the exact point of the sun in which the 
neutrino is produced. If the oscillations are of a period much shorter than this uncertainty, then the above formula 
for observing the neutrino must be averaged over one or more cycles. The result would be the apparent loss of half 
of the neutrino flux. This is not yet incompatible with the solar neutrino data. This situation occurs for md 3> 2ir in 
agreement with the condition for the validity of the simplification made in Eq. (3.10). 
Let us now consider two other ranges for md. For md <C 1, we have 



NCO 



,1/2 



/>oo 

2 d 2 / dp exp(-p 2 
Jo 



1 - - cos ( 2 ^ p 



(3.11) 



hence the conclusion that for very stringent localizations (e.g. d ~ 1/Mw) the oscillating term is negligible, as 
anticipated in the previous Section. 

When md ~ 1 (d ~ neutrino Compton wavelength), we must resort to a numerical calculation. However, in this 
case, we can restrict NCO and find its value for large time. To this aim, we set md = 1 and perform the following 
change of variable 



y = 2mt (p 2 + 1) 



1/2 



in the time dependent part of NCO, explicitly 

4e 



NCO =1 



r l/2 



dy 



2m( 



1 



Am 2 t 2 



1/2 



where 



1 



1-1/2 



dp 



1 



exp(— y 2 /Am 2 t 2 ) cost/ , 



exp(- /9 2 ) . 



With simple algebraic operations it is easy to show that NCO is limited by 



G 



T — < NCO < 1+ — . 

mt mt 

Thus, in the limit of t — + oo, we obtain 

lim NCO = 1 ~ 0.48 . 

t — >oo 

In this case the expectation value for the negative chirality average is 

(-> '-=3? 0.76 [md=l]. 

To the extent that this limit corresponds to the peak passage of the neutrino flux through the earth, roughly 25% of 
the neutrino flux would be invisible. 

The results found in this Section are shown in Fig. [I|, where NCO (and consequently (— )) is plotted for the various 
cases of md as a function of time. 
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FIG. 1. The term NCO as a function of time for md 1 (dashed), md = 1 (solid) and md -C 1 (dotted, practically zero) 
The right-hand scale represents (— ), the Average Probability of Negative Chirality (APNC), as a function of time. 



IV. CONCLUSIONS 

In this work we have suggested that a spin i particle produced in a localized condition is subject to chiral oscillations 
reminiscent of zitterbewegung. Normally this effect is of little or no interest, but for neutrinos whose absorption 
involves the V-A charged current the effect may be physically significant. We have suggested that this phenomenon 
may well explain the missing solar neutrino problem. Specifically the effect is negligible for md <g; 1. It predicts ~ 25% 
"loss" for md ~ 1, and after averaging over multiple oscillations ~ 50% "loss" for md ^> 1. With the identification of 
d of the previous Section, the critical neutrino mass m = l/d becomes 3 x 10 -15 eV. 

Only one of our basic results coincides with more traditional oscillation models: the neutrino must have a mass 
for this effect to exist. Apart from this the differences are notable. Our oscillations have periods determined, for 
large md, by the particles Compton wavelength and ignores or is in addition to any neutrino mixing. There may be 
objections to the apparent non Lorentz invariance of our final formulas, but we have already noted that the initial 
localized form is clearly frame dependent. We recall that any result specified in a particular frame is by definition 
Lorentz invariant (and can always be expressed in a manifest invariant form). Explicitly, we have chosen an inertial 
Lorentz frame centered on the sun for the initial gaussian neutrino wave function. 

We do not intend in this paper to confront all neutrino oscillation data. However, since the oscillation period is 
practically inversely proportional to the mass (see Fig. Q), we expect that the period of the muon neutrino be much 
smaller than that of the electron neutrino. Thus the atmospheric neutrino data is not incompatible with neutrino 
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masses such that the atmospheric distance corresponds to one or more oscillations for the muon neutrino while far 
too short for any significant electron neutrino oscillation. This argument could be used to determine an upper limit 
to the latter mass and the order of magnitude of the muon neutrino mass, once the localization distance d atm is fixed 
(no longer related of course to the sun's core). Rather than speculate further on this, we wish to emphasize that 
the choice of an initial gaussian probability distribution was justified by its integrability, not on the basis of physical 
arguments and should therefore be treated with caution. On the other hand the considerations in this work invite a 
serious study of the appropriate creation Hamiltonian for weak interactions so as to permit the derivation of the true 
energy eigenvalue solutions. After which a quantitative calculation, even if only numerical, can be performed. 

We conclude with a short list of the questions which this work has stimulated. First, what is the creation Hamiltonian 
of the weak interactions? How does this Hamiltonian produce a well known but often ignored fact, common to all 
created particles: the localization of the wave function? Furthermore, is it correct after the first instance to use the 
free Hamiltonian as the time evolution operator (as assumed in this work)? How are the negative frequency parts of 
the wave function to be interpreted if the reply to the previous question is affirmative? Finally, can all the known 
oscillation data be explained within the context of chiral oscillations? 
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